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Abstract 

We have modeled, tested, and installed a large, cerium-activated Gd2SiOs crystal scintillator for use as a detec- 
tor of gamma rays up to 600 MeV. We present the measured detector response to two types of incident photons: 
nearly monochromatic photons up to 40 MeV, and photons from a continuous Compton backscattering spectrum up 
to 600 MeV. Our GEANT4 simulations reproduce the measured spectra well. 
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1. Introduction 

Cerium-activated gadolinium oxyorthosilicate 
(Gd 2 Si0 5 :Ce, or GSO) HI] is a scintillator with a high 
light output and fast decay time relative to many other 
commonly used crystals, e.g. Bi4Ge30i2 (BGO). These 
properties, along with the crystal's non-hygroscopic 
nature, its relative ease of growth, and its radiation 
hardness, have made it a popular choice for a number 
of detection applications. Its most high-profile use is in 
positron emission tomography [2], where detectors are 
optimized for 511-keV photons, but GSO scintillators 
have also been used to detect protons Q3£] , charged 
leptons, andpions 0]. 

In 2009, a GSO crystal with Ce:0.5 mol% doping, 
grown by Hitachi Chemical and read out with a photo- 
multiplier tube (PMT), was adopted as a gamma detec- 
tor for the upgraded Compton polarimeter JsJ] in Hall 
A |0] of the Thomas Jefferson National Accelerator Fa- 
cility (Jefferson Lab) J3]. This device exploits Comp- 
ton scattering to make a continuous measurement of 
the longitudinal electron-beam polarization, a vital pa- 
rameter for a significant portion of Hall As experimen- 
tal program. Integration of the energy that backscat- 
tered photons deposit in this crystal allows such a mea- 
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surement to be made with accuracy better than 1 % J5t] . 
Depending on the beam energy chosen for the experi- 
ment (1-6 GeV), the maximum energy of a Compton- 
backscattered photon may range from 34 MeV to 1 GeV. 

Uncertainties in the detector response to incident 
gammas are important potential sources of systematic 
error, especially if the response is non-linear. To re- 
duce these uncertainties, a model was developed us- 
ing the GEANT4 simulation toolkit |Q] and was com- 
pared to calibration data, taken at two facilities, with 
incident gammas from 20 to 600 MeV. Section [2] de- 
scribes the Compton polarimeter and the fundamental 
simulation method as applied to single-arm Compton 
photon data at Jefferson Lab. Section [3] discusses fur- 
ther tests at Jefferson Lab using photons tagged by co- 
incident Compton-scattered electrons. In Section[4] we 
give the results of tests at the High-Intensity y Source 
(HlyS) |@] with nearly monoenergetic photon beams in 
the 20 - 40-MeV range. 

2. Simulation Method and Test with Single Photons 
at Jefferson Lab 

In the Compton polarimeter in the Hall A beamline 
at Jefferson Lab, the polarized electron beam is routed 
through a four-dipole magnetic chicane. In the cen- 
ter of the chicane, the beam passes through circularly 
polarized laser light in a high-finesse Fabry-Perot cav- 
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ity. Compton-backscattered photons pass through the 
third dipole undeflected; the GSO calorimeter is located 
on this direct path. A silicon-microstrip detector (Sec- 
tionO, located above the third dipole, detects Compton- 
scattered electrons, which are deflected through a larger 
angle than the unscattered majority of the beam. The 
Fabry-Perot cavity is periodically taken out of reso- 
nance, eliminating Compton-scattering events in the 
cavity and allowing a direct background measurement. 
During typical running, the dipole fields are controlled 
via feedback from a beam-position monitor in the chi- 
cane, so as to maintain a stable beam position. 

Light from the GSO calorimeter is collected in a 2- 
inch, 12-stage PMT, with a base customized for max- 
imum linearity of response. The data-acquisition sys- 
tem is based on a 200-MHz, 12-bit Struck SIS3320 flash 
analog-to-digital converter, modified to integrate the in- 
put signal over an externally timed window. In addi- 
tion to performing this onboard integration, the SIS3320 
card records all samples from the window in one of two 
internal buffers. Timestamps from photon- or electron- 
detector triggers, recorded in a CAEN V830 latching 
scaler, allow the retention of some information about 
individual pulses. For a prescaled sample of pulses, the 
numerical sum of samples from the programmable read- 
out window is recorded in the datastream; the energy of 
the incident photon can be retrieved from this pulse in- 
tegral. For a smaller number of pulses, all of the sam- 
ples for the readout window are written to disk, allow- 
ing pulse-shape analysis J^t]. It is also possible to trig- 
ger the system at regular intervals that are uncorrelated 
with photon pulses, allowing the study of pileup events. 
These triggers are generated either in software or with a 
remotely programmable function generator. 

The GEANT4 Monte Carlo (MC) simulation of the 
GSO crystal response begins with the generation of a 
beam of simulated photons to match the experimental 
beam. To reproduce a Compton-backscattered beam, 
simulated photons of various energies are generated 
with probabilities weighted by the Compton scattering 
cross-section for the specific kinematics (i.e. initial elec- 
tron and laser-photon energies) of the planned experi- 
ment. The simulated photons are then allowed to inter- 
act with items in the beamline; in the standard Hall A 
installation, this includes a 1-mm-thick, 4-cm-diameter 
lead filter and a 5-cm-thick, 8-cm-diameter lead colli- 
mator with an interchangeable aperture up to 2 cm in di- 
ameter. The final item in the beamline is the GSO crys- 
tal, a cylinder 6 cm in diameter and 15 cm (10.9 radia- 
tion lengths) long; TableQ]lists the other crystal proper- 
ties used in the MC. Figure [1] histograms the locations 
where photons in the simulation first interact with mat- 
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Figure 1 : Histogram of positions of the first interaction point of inci- 
dent photons in the Hall A Compton gamma beamline, as determined 
by the GEANT4 MC. The z axis is the central axis of the incident 
Compton-scattered photons; the y axis is vertical. In this figure, the 
central axis of the calorimeter and the 2-cm aperture of the collimator 
is offset by 0.5 cm from the gamma beam, reflecting conditions during 
one run period. 



Property 



Density 
Radiation length 
Attenuation length 
Birks' constant 



Value 

6.71 g/cm 3 
1.38 cm 
340 cm 
5.25 Lim/MeV 



Reference 



[12] 
□J 

estimated from [3] 
Li3J 



Table 1: GSO:Ce (0.5 mol%) properties used in simulation 



ter. GEANT4 modeling of the gamma shower relies 
on the cross-section cr for gamma conversion into an 
(e + , e~) pair; for gamma energies between 1.5 MeV and 
100 GeV, the parameterization is accurate to within 5%, 
with a mean accuracy of 2.2% d. Electromagnetic in- 
teractions are modeled based on the Livermore physics 
list; hadronic interactions are based on the QGSP_BIC 
physics list. The combined non-linearity of the pho- 
tomultiplier tube and front-end electronics was deter- 
mined empirically and the resulting functional 
form is an input to the simulation. 

An optical extension to this basic simulation fol- 
lows the path of each scintillation photon produced in 
the electromagnetic shower from an incident Compton- 
scattered photon. This package approximates the pol- 
ished GSO surface as perfectly smooth, and includes 
the aluminum-foil detector wrapping (with a modeled 
reflectivity of 0.9) and the efficiency of the PMT photo- 
cathode. The output of the MC is the number of photo- 
electrons produced in the simulated photocathode. The 
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Figure 2: A Compton spectrum measured in Hall A, with the ADC re- 
sponse given in arbitrary units. The Compton edge in this configura- 
tion was at 204 MeV. A non-optical MC with 2.3% Gaussian smearing 
is fit to the experimental data with only two free parameters: a hori- 
zontal scale factor and a vertical scale factor. The discrepancy at low 
photon energies is an artifact of an anomaly in the trigger threshold; 
this region is excluded from the fit. 



resulting spectrum shows non-Gaussian smearing due 
to optical effects such as shower leakage from the crys- 
tal and the dependence of photon-collection efficiency 
on the initial interaction position. However, the com- 
putational cost of this package is prohibitive for some 
simulation tasks. 

Spectra simulated under Jefferson-Lab conditions 
must be modified to take into account pileup, in which 
two or more pulses arrive and are integrated during the 
readout window for a single trigger. This correction is 
determined experimentally based on periodically trig- 
gered snapshots of the readout window, which give the 
random-event rate and the spectrum of energy deposited 
by random pulses during such a window. For spec- 
tra from electron-photon coincidence data (Section [3J, 
it is sufficient to use this spectrum to add random, 
pileup pulse integrals to those of simulated Compton- 
scattered photons. For photon-arm singles data, how- 
ever, the correction is complicated by the possibility 
of pileup between two Compton events; a naive pileup 
correction would count such a pair twice. This double- 
counting effect is canceled by a careful construction of 
the empirical pileup spectrum that is added to the MC: 
half from dedicated background measurements and half 
from Compton-scattering data 1 1 4TI . For simulations of 
the primary running mode, in which signals are inte- 
grated over a long, untriggered window, no pileup cor- 
rection is necessary. 

A 2.3% Gaussian smearing function must be folded 
into the basic, non-optical MC for a satisfactory match 
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Figure 3: Fits of the non-optical MC, including background, to 
tagged photon spectra fr om |(a)| strip 10 (125.1 - 127.0-MeV pho- 
tons, x 2 /dof = 2.13) and|(b)|strip 21 (158.5 - 160.4-MeV photons, 
X 2 /dof = 2.33). 



to the experimental data, as shown in Fig. [2] By 
contrast, the optical MC requires only 1.5% Gaussian 
smearing; that is, a 1.5% smearing is not accounted for 
by modeled optical effects. At high data rates, pileup 
corrections are necessary for both models. 

3. Tests with Tagged Photons at Jefferson Lab 

Further tests of the GSO calorimeter in Hall A were 
performed with photons tagged by the electron detector 
located above the third dipole of the Compton polarime- 
ter. Each of this detector's four planes consists of 192 
horizontal strips with a 240- [im pitch. When a scattered 
electron strikes the detector, the vertical position of the 
strip gives its deflection angle in the third dipole and 
therefore its momentum, which in turn determines the 
energy of the associated scattered photon. The electron 
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Figure 4: Measured spectrum of pileup pulses in tagged-photon data, 
used to correct MC spectra. 



detector can thus tag the energies of coincident photons 
seen in the GSO crystal. For these data, taken with an 
initial electron-beam energy of 3.4 GeV, only the bot- 
tom 37 strips saw Compton-scattered electrons. 

Figure [3] shows measured and modeled spectra of 
photons tagged by two typical strips. GSO readout was 
triggered by the electron detector. Due to computational 
constraints, the MC is non-optical, but it does include 
simulated background events and pileup effects; Fig. [4] 
shows the measured pileup spectrum used for the MC. 
Each fit has three free parameters: a vertical scale factor, 
a horizontal scale factor, and a factor setting the amount 
of included background. The latter varies from strip to 
strip because electron-detector noise, which is slightly 
different for each strip, is a significant source of back- 
ground. MC fits for all 37 strips are given a common 
horizontal offset and 5% Gaussian smearing; the larger 
smearing requirement is likely due to poorly understood 
effects in the electron detector. 



4. Tests with 20 - 40-MeV Photons at HlyS 

HlyS produces a 2- to 65-MeV gamma beam via 
Compton scattering. At HFyS, electron bunches that 
circulate in a storage ring are sent through magnetic un- 
dulators, causing the electrons to produce free-electron 
laser (FEL) light, which is stored in an optical cavity. 
The FEL photons are then allowed to collide with elec- 
trons in the storage ring; the Compton-backscattered 
photons travel through a manually adjustable collimator. 
The result is a nearly monoenergetic gamma beam in 
the experimental hall, approximately 60 m downstream 
of the Compton interaction point. 
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Figure 5: Simulated spectra fit with two free parameters (horizontal 
and vertical scale factors) to experimental data from the HlyS facil- 
ity, for incident photons at |(a)| 20 MeV (x 2 /dof = 3.26, plotted with 
the ratio of measured counts to simulated counts (dashed line)), |(b)| 
22 MeV (x 2 /dof = 10.97), and[p]25, 30, and 40 MeV, with a com- 
bined fit for all three energies (x'/dof = 10.20). Gain shifts prevent 
combining the fits to the 20- and 22-MeV spectra. 
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Figure 6: The measured Compton asymmetry plotted as a function of 
ADC response to deposited photon energy. The experimental data are 
compared to GEANT4 MC data with no adjustable parameters. The 
Compton edge in this configuration was at 204 MeV. 



A series of HlyS-beam test runs at five different pho- 
ton energies (20, 22, 25, 30, and 40 MeV) was taken us- 
ing the GSO calorimeter; gain shifts after the 20-MeV 
and 22-MeV data points prevent the analysis of all five 
energies together. Figure [5] shows the measured spec- 
tra at these energies, as well as the ratio of measured to 
modeled counts for 20-MeV photons. An optical MC of 
the HlyS configuration, including non-linearity, was fit 
to each spectrum with two free parameters: horizontal 
and vertical scale factors. Background and pileup ef- 
fects were negligible at all energies except 22 MeV, and 
were not included in the MC. 

Because the energy spread of the HlyS beam is well 
known, fits to these data constrain the minimum amount 
of additional smearing that must be included to bring the 
MC into agreement with the data. The corresponding 
1.5% smearing is included in the simulations shown in 
Fig. 13 



5. Conclusion 

We have demonstrated agreement at the few-percent 
level between the MC of this GSO:Ce crystal and pho- 
ton spectra measured over a range of rates and en- 
ergies. The ultimate purpose of the MC, however, 
was to unfold the detector response from Compton- 
photon polarimetry data in Hall A. The asymmetry in 
Compton-scattering cross-sections for opposite photon- 
and-electron spin configurations is proportional to the 
polarizations of both incident beams and to the analyz- 
ing power of the apparatus |5J]; the MC is a vital in- 
put to the latter. Figure [6] shows the measured Compton 



asymmetry as a function of the energy deposited in the 
GSO by Compton-scattered photons. The non-optical 
MC result is plotted on the same axes with no adjustable 
parameters. The horizontal scale was determined from 
a fit to the full photon energy spectrum (Fig. [2}; the 
vertical scale is set by the electron beam polarization, 
which was confirmed by M0ller polarimetry measure- 
ments [ 15], and by the incident photon polarization. 

With the aid of this MC, the analyzing power for the 
GSO calorimeter is known to within 0.33%; this un- 
certainty is dominated by the 0.3% error on the detec- 
tor non-linearity measurement Q]- Existing GEANT4 
physics modules allow simulation of the GSO:Ce re- 
sponse to gammas in the 20 - 600-MeV energy range 
with satisfactory precision. 
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